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We calculate the AC linear response of a superconductor in a nonequilibrium electronic state. 
The nonequilibrium state is produced by injecting quasiparticles into the superconductor from nor- 
mal leads through asymmetric tunnel contacts. The dissipative part of the response is drastically 
increased by the injected quasiparticles and is proportional to their total number regardless of the 
imbalance between the numbers of electron-like and hole-like excitations. 

PACS numbers: 74.25.N-,74.40.Gh 



I. INTRODUCTION 

Over the years, measurement of the impedance was 
used extensively in the studies of spectra of elementary 
charge carriers in normal metals and superconductors. 
In the case of normal-metal single crystals, measure- 
ments of surface impedance in conditions of cyclotron 
resonance were instrumental in the reconstruction of the 
Fermi surface geometryii In superconductors, the tem- 
perature dependence of the impedance allowed one to in- 
vestigate the appearance of the BCS gap in the spectrum 
of quasiparticlesj^ Recent interest to the AC response of 
superconductors is driven by its application for studying 
new materials^i^ and by the use of well-studied supercon- 
ductors as material for high-quality resonators. Depend- 
ing on the problem, it is important to realize the lowest- 
possible dissipation rate of a resonator, or its well-defined 
response to a perturbation. The former goal is central for 
superconducting qubit physics^ and quantum optics^^ 
where one is interested in achieving the longest-possible 
coherence times. The latter one is important, e.g., for the 
use of superconductors in microwave kinetic inductance 
detectors j^ 

To interpret the measured impedance, one has to com- 
pare it with theoretical predictions. However up to now, 
no microscopic calculations of impedance of superconduc- 
tors under nonequilibrium conditions were available. The 
frequency and temperature dependence of the complex 
conductivity of a superconductor within the BCS theory 
was first evaluated in the seminal paper of Mattis and 
Bardeeufi^ where the limit of short electron mean free 
path (the "dirty superconductor" limit) was considered. 
The non-local conductivity of a clean superconductor was 
derived by Abrikosov, Gorkov and Khalatnikoviii Later 
works of Naroi^ were aimed at bridging the two limit- 
ing cases. The theory developed in Refs. [l0l - [l2l addresses 
superconductors at thermal equilibrium. If desired, this 
condition is fairly easy to achieve at not-too-low tem- 
peratures. Indeed, even the early measurements have 
demonstrated the hallmarks of the BCS behavior of the 
dissipative part of the impedance, including its thermal- 
activation behavior at sub-gap frequencies J^ At the same 
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FIG. 1: N-I-S-I-N structure: the superconductor (S) is con- 
nected to two normal leads (N) - maintained at different volt- 
ages (bias V = Vr — Vl) - via tunnel junctions characterized 
by tunneling rates Fr, Tl- The superconductor is also sub- 
ject to a weak electromagnetic field oscillating with frequency 



time, experiments involving impedance measurements in 
non-equilibrium conditions relied on heuristic extension 
of Mattis-Bardeen formula, amounting to the replace- 
ment of the equilibrium quasiparticle distribution func- 
tion in the formula by a non-equilibrium one.^^^ 

The goal of this paper is to microscopically evaluate the 
linear AC conductivity of a superconductor in a concrete 
setup allowing a controlled perturbation of the quasipar- 
ticle distribution function. To this end, we analyze the 
steady state and linear AC conductivity of a supercon- 
ductor brought out of equilibrium by electron tunneling 
through two junctions connecting the superconductor to 
normal leads (N-I-S-I-N structure - see Fig.[T]); the struc- 
ture is biased by a constant voltage V. 

We find the steady-state electron distribution at finite 
temperature T and voltage V, assuming these two scales 
small compared to the quasiparticle energy gap A. In the 
case of unequal conductances of the two tunnel junctions, 
charge imbalance is created along with a perturbation of 
the energy distribution of quasiparticles. We evaluate 



the complex AC conductivity of the superconductor in 
this non-equihbrium state and show that only the en- 
ergy mode of the quasiparticle distribution enters in the 
proper generalization of the Mattis-Bardeen formula. At 
eV, ksT <ti A we cast the result for conductivity in terms 
of T and quasiparticle density nqp. That form extrapo- 
lates between the equilibrium result (where Uqp is a func- 
tion of T only) and the non-equilibrium one, where Uqp 
is a function of V and T, with arbitrary ratio eV/ksT. 

In the next Section, we formulate the problem in terms 
of matrix Green functions. The steady state of the elec- 
trons in the superconductor formed in the presence of 
finite bias applied to the N-I-S-I-N structure is found in 
Section IIIH where we also establish the correspondence 
between the descriptions in terms of the matrix distribu- 
tion function F{e) and the scalar distribution function 
of quasiparticlcs /e (the latter was used in the original 
Mattis-Bardeen theory for the AC conductivity at equi- 
librium). Using the description in terms of /^, we in- 
clude the electron-phonon interaction into the consider- 
ation of quasiparticle kinetics. The AC conductivity at 
low temperature and bias (but at arbitrary eV/kBT) is 
analyzed in Section IIVI Throughout the paper, we use 
units h = kB = i- 



II. ELECTRON DYNAMICS IN A 

SUPERCONDUCTOR SUBJECT TO DC BIAS 

AND WEAK AC FIELD 

We consider a diffusive superconductor connected to 
two normal leads, left (L) and right (i?), via tunnel bar- 
riers and exposed to an external, time-dependent electric 
field. The system properties can be described in terms of 
disorder-averaged matrix Green's functions for the super- 
conductor, ^s(k, t, t'), and for the electrodes, Gi{pi, t, t'), 
i = L,R. Each of the matrices Q has the form, in Keldysh 
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The elements of this matrix are 2x2 matrices in Nambu 
particle-hole space. The superconductor Green's func- 
tion obeys the Dyson equation 
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Here A(t) is the vector potential, which is related to the 
electric field via E — —dA/dt. The coefficient l/r is the 
impurity scattering rate, and vs is the density of states 
at the Fermi level in the superconductor. The matrix 



element Ti for tunneling into lead i determines the di- 
mensionless conductance gi = ^ir'^vsViTf of junction i, 
where vi is the density of states in the lead. The assump- 
tion 



gL +gR < gs, 



(3) 



where gs is the normal-state conductance of the super- 
conductor, justifies the use of the tunneling Hamiltonian 
from which the last term in Eq. ([2]) is derived. The same 
assumption enables us to neglect small spatial variations 
of the order parameter, which we take to be uniform, real, 
and time- independent. The last two conditions amount 
to a choice of gauge. 

Since we are interested in the linear response to the 
external field, we first consider the system in the absence 
of field, but with the leads at different potentials. This 
is the focus of the next section. 



III. NON-EQUILIBRIUM STEADY STATE 

When there is no external field (A — 0), the Green's 
functions are isotropic in momentum space and depend 
on the the difference t—t'. Then the semiclassical Green's 
functions 
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where ^p = p^ /2m — Ep, depend only on the time differ- 
ence, and not on the momentum direction n. Therefore, 
we omit the argument n for the rest of this section. In 
the normal leads, the elements of d are (in the frequency 
domain) 



Ct,- 



Gr, 



^i — ^i '^i ~ "^i^i 1 



(5) 
(6) 



with 



n^(e) 







, , Tr\ \ , n(£) = tanh, 
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(7) 
The potentials Vi in the leads are measured from the su- 
perconductor chemical potential. Using Eq. (O, Eq. (|6]) 
can be rewritten as 



where 



Gf = 2^0 f,+2nii. 






(8) 
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The superconductor Green's function Gs is determined 
by Eq. ©; taking the difference between Eq. ([2]) and its 
conjugate, and using Eqs. (|4|), we find 

[sT, + iAfy, Gs] = -I [Tl Gl + Tr Gr, Gs] , (10) 



where Fj = nVi Tf are the tunnehng rates and Vi are the 
densities of states in the leads i = L,R. This equation 
must be supplemented by the normalization condition,— 
which in the frequency domain takes the form 



If |e| < A, the sum ^^ + ^"^ is of the order of A, i. e. 
much larger than Tl and F/j. However, if \e\ > A this 
quantity is proportional to these small parameters and 
can be approximated as 



Gs{e)Gsis)^i. (11) 

Isolating the R{A) component of Eq. ([TU| gives 

where 



^mA) ^ (^ ^ .p^ ^ .p^) ^^ ^ .^ -^ 



In view of the R{A) component of the normalization con- 
dition dni), (^5*^')^ = i, the solution to Eq. ^ can 
be written as 
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where 



ili(^)=±[(s±^TL±^Tnr~A'] 



1/2 



(14) 



(15) 



In the limit F, -^ 0, Eq. (|T4|) reduces to the well-known 
semiclassical BCS expression.— When Fi 7^ the exci- 
tations in the superconductor have finite lifetime due to 
tunneling into the normal leads and this causes broad- 
ening of the density of states.— In what follows we as- 
sume that the broadening is much smaller than the gap, 
Fl + Ffl <C A. Since the dimensionless conductance, as 
defined in the text before Eq. (O, equals gt = SnTi/S, 
with 6 being the level spacing in the superconductor, we 
can express this condition as 



gL + gR < A/5. 



(16) 



Wc now consider the Kcldysh component of Eq. ([TU)) . 
It may be written in the form 



^« G§ -GlH^^ -iTl [Gf Gj 
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Its left-hand side can be simplified as follows: first, we 
use Eq. (|14p to rewrite H^'^^^ in terms of Gg ; then we 
replace the advanced Green's function with the retarded 
one employing the orthogonality condition [the Keldysh 
component of Eq. pT|)] 



GR /^K I /^K /^A 
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Multiplying the resulting equation by Gg from the left 
and using again {GgY = 1, we arrive at 
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where 



e^(^^ = ±[(£±io+)-A2] 



1 

'¥ 
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Consistently with this approximation, the Green's func- 
(13) tions Gg in the right-hand side of Eq. (fT7|) should be 



taken at zero order in Fj's. Using Eqs. (|5]) and ([H]), after 
some algebra we obtain 



Lfg — 



1 1 



n°{ef^ +iAfy)+ no- 



where (a — 0, 1) 



FL< + Ffl,n^ 
Tl+Tr 



(21) 



(22) 



with n1 j^ of Eq, 



dH). We stress that Eqs. P^ - (PT|) 
are not restricted to junctions with conductances smaller 
than the conductance quantum, the conditions of appli- 
cability being given by Eqs. ^ and (fT6| . 

The matrix G^ depends on the voltages Vl, Vr via 
Eqs. dl]) and p2|). However, under steady-state condi- 
tions these two quantities are not independent and, as we 
discuss in the next section, the state of the superconduc- 
tor is fully determined by their difference V — Vr — Vl- 



A. Current and potentials 

The results of the previous sections rest on the steady- 
state assumption. It means that the total current fiowing 
out of the superconductor must vanish. This require- 
ment, as we show below, defines the division of the ap- 
plied bias between the two tunnel junctions connecting 
the superconductor to the normal leads. 

The current leaving through the left tunnel contact is 
given by 



^L = ^ ei^sFi 



deTr 



1 nR nK I /^K AA 

T^l U^ Ug + U^ Ug 



G" R /^K /^K /^A 
S ^L ~ ^S ^L 



(23) 



The current Ir into the right contact is found by replac- 
ing L ^f R. Therefore the total outgoing current is 



1 
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Il+Ir^evs / d£\ -^ - -J ] —^ [Vl n 



llL 



Fflnifl) 
(24) 



and must be zero under steady-state conditions, as dis- 
cussed above. [Wc used Eqs. ([5|), p^ . and (PT|) in deriva- 
tion of Eq. ()24|) .] In other words, the requirement 



Il + Ir = 



(25) 



determines the relationship between the voltages Vl, Vr, 
that enter as parameters in Eq. ([33]). 

Let us consider explicitly the regime of low tempera- 
tures and voltages so that maxjeVi, cVr, T} <g; A. Then 
the following approximate expression for the functions 



'L(fl) 



are valid at energies |e| > A: 



nl{e) ~ -2 sinh(eVi/T) exp(-|e|/T) . (26) 

Substituting into Eqs. p4 )) - ((25)) we obtain the condition 

Tl smh{eVL/T) + Tr smh{eVR/T) = 0. (27) 

Introducing the bias V — Vr — Vl, we find from this 
equation that 



smh{eVL/T) = - 



Tr smh{eV/T) 



V^l +rl + 2r lTr cosh{eV/T) 



(28) 

Thanks to this relation, we can reexpress voltages Vl , Vr 
in terms of the bias V, and the non-equilibrium steady 
state in the limit eV, T <C A is fully determined by these 
quantities. 



B. Quasiparticle distribution function 

In the preceding sections we have found an explicit ex- 
pression for the non-equilibrium Keldysh Green's func- 
tion for the biased superconductor, Eq. ((2T|) . In this 
section we derive the relationship between the Keldysh 
Green's function and the distribution function / for the 
Bogoliubov quasiparticles. 

Our starting point is the definition of Gg (in the time 
domain) in terms of creation and annihilation operators 

G§{t,t') (29) 



-nvs 



Z^ 



The sum over single-particle states — ^^ with energy 
S,m is equivalent to the integration over their energy ^p 
in Eq. ((4]). Next, we perform the Bogoliubov transfor- 
mation from electron operators Cm-\, c^i to quasiparticle 
operators am, l3m'- 

with amplitudes Um, Vm given by 



^m^n 



1 A 
' ~ ~2^ 

Zj C771 



(31) 



and em = \/^m + ^^ being the quasiparticle energy. 
Introducing the distribution function (wc assume equal 
population of the two spins) 



fiU 



\^m ^n 



(/3,l/3. 



(32) 



using Eq. pi|) . and taking the Fourier transform with 
respect to t — t' wc obtain 



Gf (e) = 2Re^=L=={(£f, + *Ar,)[l - 2/^(£)] 



2^/1 



'f^ie)} 



(33) 



with the energy and charge modes of the distribution 
function f{Cm) defined as 



/^(e) = - /(Ve2-A2) + /(-V^ 



/<3(e) = - [/(v/e2-A2) - f{-Ve^~A^) 



(34) 



Hence the Keldysh Green's function can be presented in 
the form 
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where 



F = sgne[l-2/^(e)]i 



V? 



Tz +i- 



(35) 



2/Q(e). (36) 



This shows that the trace of F is directly related to the 
energy mode / , while the traceless part possesses an 
additional density of states factor compared to the charge 
mode f^ of the quasiparticle distribution function. 
Comparison of Eq. (p6| with the solution of the kinetic 



equation for Gg, Eq. (|2T|) . enables us to find the explicit 
expressions for the distribution function modes 

1 






[1 -sgnen"(e)]. 



(37) 



and the voltages Vl, Vr are related by Eq. ([25|) . At low 
temperature and bias T, eV <C A, at leading order in 
e~^/'^ we obtain, using Eq. p8|) . 



/^(e)«<f>(ey/T)e--/^, 



(38) 



with 
<^{eV/T) = 



Tl + rl + 2rLTR cosh{eV/T) , 



(39) 



while at this order f^ vanishes due to Eq. ([37)1 . 

One may also obtain the distribution function /c by 
solving the kinetic equation for quasiparticles. The lat- 
ter method is convenient to assess the effects of relax- 
ation. In the next section we consider in more detail the 
electron-phonon interaction. 



C. Kinetic equation for quasiparticles 

In the kinetic equation approach, the effects of inter- 
action are described by adding to the rate of change of 
the occupation numbers due to tunneling an appropri- 
ate colhsion integral /. In the case of electron-phonon 
interaction, one can derive a set of coupled kinetic equa- 
tions for the two modes of the quasiparticlc distribution 
function, f^ and f^: 






dV 



TlIol + Tfl/ofl - r/^ (40) 

/f{/^,/'=,iV}+/f{/^,/«,iV}, 



Ve^ - A2 



[rL/iL + rfl/iflj-r/Q (4i) 



+ iQ{f^j^,N}+iQ{f^jQ,N}, 



where T = Tl + r^, /o, = (1 - nO)/2, and fu = nl/2. 
Neglecting the collision integrals, in the steady-state 
df/dt = we recover immediately the result Eq. ([37| . 
The subscripts r and s arc used to distinguish, in the 
electron-phonon collision integrals, quasiparticlc recom- 
bination and scattering processes, respectively, and we 
take the phonon distribution function N to be the ther- 
mal equilibrium one: 



N, 



1 



.J/T _ I 



(42) 



We characterize the electron-phonon collision rate 
by the scattering rate, l/rph, for normal electrons off 
phonons at the critical temperature Tcr^ We consider 
the limit of low temperature and bias, in which case 
the electron-phonon interaction has a small effect on the 
quasiparticlc distribution as long as tunneling is the main 
scattering process. Considering relaxation rate for an ar- 
bitrary small perturbatioi>i^ of the quasiparticle distribu- 
tion function, one might conclude that the correspond- 
ing condition for that is {l/Tph){T/Ay/'^ < T. (The 
small factor (T/A)^/^ here suppresses the low tempera- 
ture scattering rate for quasiparticles compared to that of 
normal electrons at T^.) However, for the specific case of 
out-of-cquilibrium distribution created by tunneling the 
condition for the rate l/Tph is even softer, sec Eq. ((50l) 
below. 

Evaluating the correction to the charge mode f^ due 
to the weak electron-phonon relaxation at low tempera- 
tures, we find that such correction to Eq. ([37| vanishes 
at leading order in e~^l'^ due to the relation Eq. (P7)) . 

We now consider the energy mode. The collision inte- 
grals ifg contain terms quadratic in /'5, which we can 
neglect since /*3 vanishes at leading order, as discussed 
above. The non-vanishing terms are^S 
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for recombination processes and 
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(44) 



for scattering processes. 

We solve the kinetic equation (^0)) by iterations, writ- 
ing /^ in the form: 



r-/ 



E(0) 



1 



Tp/iF 



Ell) 



(45) 



with /-^(o) given in Eq. ([55]). At leading order in e~^/'^ , 
only the scattering collision integral if is present. It 
is satisfied by any Boltzmann distribution function, like 
the one in Eq. (p8)) . Indeed, let us consider the terms in 
square brackets in the second line of Eq. (|44| ; in the low 
temperature and bias regime they are approximately 



$e~'''/^7Ve_e' - $e~''/^(7Ve_^' + 1) = , 



(46) 



in agreement with the detailed balance in the absence of 
recombination. Similarly, the terms in square brackets 
in the last line of Eq. ((44]) add up to zero. Therefore, at 
order e~ ' "^ there arc no corrections to the distribution 
function. 

To calculate the correction due to recombination, we 
note that e + e' > 2 A, so that in the last term in square 
brackets in Eq. (1^5]) we can neglect N^+^i in comparison 
to unity. The square brackets then becomes approxi- 
mately 



-{e+e')/T 



1_$2 



(47) 



i.e., they are of order e 2A/t^ p^j. ^^^^ correction we find 



1 r°° 



x(e + e')^(l + ^ 



A Ve'2 - A2 

-e'/T 



(48) 



After the substitution e' = A + Tx, the integral can be 
evaluated to give, at leading order 



/ 



Eil) 



g-s/T ^^ _ ^2^ 



-A/T 




[Note that at zero bias $(0) = 1 and the correction van- 
ishes.] From this equation we can estimate the ratio be- 
tween the leading term and the correction in Eq. (j45| 
and find that the iterative solution is applicable if 




A/T 



(50) 



Under this condition, the electron-phonon relaxation 
leads only to a small modification of the distribution 
function given by Eqs. (|45|) and (|49)) . Albeit the cor- 
rection Eq. P^ at low energies {e « A) scales with tem- 
perature predominantly as e~^ '^, it may exceed the 
corresponding term of expansion of Eq. (|37p in e^ '^ . 
That happens at large values of l/rp^F, which still may 
satisfy Eq. ([SO)) . 

This concludes our analysis of the non-equilibrium 
steady state of a superconductor for the specific prob- 
lem of a finite-bias tunnel injection of quasiparticles. In 
the next section we consider the linear response of an out- 
of-equilibrium superconductor to an external AC field. 



IV. AC CONDUCTIVITY OF N-I-S-I-N 
STRUCTURE 

In this section we study the linear response of the su- 
perconductor to an external AC electric field oscillating 
at frequency uj. Within linear response, the Green's func- 
tion Gs is the sum of the zero-approximation isotropic 
part Go and a small correction Gi, linear in A, which is 
anisotropic in momentum space 

G5(n)-Go + Gi(n), n = p/|p|. (51) 

The correction Gi(n) determines the current density via 



The zero-order in Ti approximation is valid when T^ + 
F/j <C min{l/T, w. A}. In particular, the condition 
Fl + F/f <C 1/r implies that the tunneling has a neg- 
ligible influence on the quasiparticle states; in fact, the 
only role of the leads is to generate a non-equilibrium 
quasiparticle population, and the calculation of the AC 
response docs not depend on the specific way in which a 
non-equilibrium population is established (as long as it 
does not substantially alter the quasiparticle states). 

Equation ([SS]) can be solved for a generic relation be- 
tween 1/r and A - see Appendix |^ Here we consider 
the dirty limit tA ^ 1. To calculate the current, we only 
need the following expression for the Keldysh component: 



Gf = ievA^T 
X [g^ f, (G^' ^G^)F^+F+ (G^ 



Gi)T,G^ 



(57) 



with F± = F{e ± to/ 2). The retarded and advanced 
Green's functions in this equation are given by Eqs. (J13p - 
(ITil) . Substituting Eq. dST]) into Eq. ^^, we arrive at 



j,^ = fT(w)E^ 



(58) 



where cr(w) is the complex conductivity. It is convenient 
to isolate its equilibrium zero-temperature kinetic part 



t7o(w) = - 



TTCTArA 



(59) 



^ei^s I deTr{vT^G^ 



K\ 



(52) 



where v ~ vpn, vp is the Fermi velocity, and angular 
brackets denote averaging over the momentum direction. 
We emphasize that the calculation of the AC response 
that we present here is not restricted to the particular 
non-equilibrium state considered in the previous section; 
rather, it is valid for a generic matrix distribution func- 
tion F , which determines the Keldysh part of the zero- 
approximation Green's function Go via Eq. p5p. 

The two terms in the sum ((5T|) obey the orthogonality 
condition 



where 



G+Gi(e,cj) + Gi(e,w)G_ = 0, 



G± = Go(e±c^/2). 



(53) 



(54) 



An equation for Gi can obtained as before by consider- 
ing the difference between Eq. ([2]) and its conjugate and 
using Eqs. ^. Assuming A{t) = A^^e"*"* we find, at 
zeroth order in the tunneling rates F^ 

-iH+Gi+iGiH^ + ^G+Gi-^GiG- ^^^^ 
It It (55) 

= ievA„(f^ G_ - G+ fz), 
where 

H±^{e±uJ|2)T^^-ib.Ty. (56) 



which represents the purely inductive response of the su- 
perconducting condensate to the external field in the ab- 
sence of quasiparticles and where ajq = 2e^vsD is the 
normal-state conductivity of the superconductor, with 
D ~ VpT/3 the diffusion constant. Using the relation- 
ship Eq. pop between matrix and quasiparticle distribu- 
tion functions, assuming < lu < 2A, and performing 
simple rearrangements, we find 



a{uj) = ao{uj) + 



de ■ 



2crjv 



£ (e + w) + A2 



A 



v/(e2-A2)[(£ + w)2-A2 



A+ui 
A 






(60) 
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de ■ 



^(e2-A2)[A2-(e-Lj)2 
e (e - w) + A2 



A + LO 



^(£2-A2)[(e-a;)2-A2] 



f^ 



Equation (|60| may be viewed as a generalization of the 
Mattis-Bardeen formula to an arbitrary quasiparticle dis- 
tribution function /j. 

We now use Eq. ([50)) to find the AC conductivity for 
the biased N-I-S-I-N structure considered in Sec. lIIIl con- 
centrating on the limit of low temperatures and voltages, 
inax{eVL, eVii,T} <g; A. In this limit, f^ is exponen- 
tially small and approximately given by Eq. psp . If in 



addition to the assumptions of low temperature and volt- 
ages made above, the frequency is also much smaller than 
the gap, cj ^ A, we find at leading order in A 



a{uj) 



X 



l-7r + 2e~^/^<P{eV/T) 



7rexp(-a;/2T)/o(w/2r) 
2iamh{u}/2T)Koiu}/2T) |, 



(61) 



where /q and Kq denote the modified Bessel functions of 
zeroth order. 

The effect of the finite bias on the AC response of the 
superconductor described by Eq. (|61]) is contained in full 
in the function $, Eq. (|5U)) . In fact, as we now show, 
this function accounts for the non-equilibrium density of 
quasiparticles. The quasiparticle density is given by [cf. 
Eq. dSlD] 



2vs jd^fiO^^^s j de 



Ve2 - A2 



/f (62) 



with /■^ defined in Eq. (|34|) and the factor 2 account- 
ing for spin. At low temperature and voltages, we can 
approximate /^ as [see Eq. (f38| ] 



/f ~ ^{eV/T)e 



-e/T 



(63) 



and evaluating the integral at leading order in A we arrive 
at 



= 4i/s$(eF/T) 



ttAT 



-A/T 



(64) 



Using this equation, Eq. (|6ip can be written as 
a{uj) = a'{u) +icr"(w), 



(65) 



cr'(w) 
a"{uj) 



2A 



(Tat — sinh —-Ka[ — 1 
u 21 



ttA 



gn- 



1 



\2TJ V 2ttTvsA' 
'\2TJ V 2ttTvsA 



We note that the above result for the AC conductivity 
in terms of the quasiparticle density preserves its form 
even if we include the correction to the distribution func- 
tion due to the electron-phonon interaction [assuming the 
condition Eq. (j50p is satisfied] . Of course, the quasipar- 
ticle density Eq. (|64p should be corrected to account for 
Eq. (|49| ^ Equation (|65|) agrees with the results obtained 
in Ref. ^ via a phenomenological generalization of the 
Mattis-Bardeen formula that includes an effective chem- 
ical potential for quasiparticles. 



V. SUMMARY AND DISCUSSION 

In this work we have derived an extension of the 
Mattis-Bardeen formula for the AC conductivity of a 



dirty superconductor to include a non-equilibrium occu- 
pation of the quasiparticle states. Equation (j60p shows 
that the charge imbalance does not affect the linear re- 
sponse to external radiation. 

As an explicit example, we considered properties of 
a superconductor connected via tunnel junctions to two 
normal leads. Application of a constant bias to the leads 
creates a steady non-equlibrium state in the supercon- 
ductor, see Sec. IIIII We have evaluated bulk AC con- 
ductivity a{ijj) of such N-I-S-TN setup for low frequency 
a; ^ A of a weak external field and low temperature 
(T < A) while keeping fixed the DC bias V. The role 
of the finite bias is to create a non-equilibrium quasipar- 
ticle density Uqp. Remarkably, it is possible to express 
the conductivity cr(cj) as a function of temperature and 
Uqp only, see Eq. (|55|) . The dependence of 17(0;) on bias 
V enters only through the 1^-dependence of Uqp. That 
simplification works even if electron-phonon relaxation is 
taken into account [if the condition Eq. ([50]) is satisfied] . 
At zero bias {V = 0) the distribution Uqp becomes the 
equilibrium distribution function, and Eq. (j65p reduces 
to the conventional Mattis-Bardeen formula. 

As a further extension one may include the electron- 
electron interaction. If sufficiently strong, it would lead 
to the replacement of temperature T in the quasiparti- 
cle distribution function with an effective temperature 
Teff. Interestingly, for a fixed non-equilibrium quasipar- 
ticle density, an increase of temperature (or effective tem- 
perature) leads to a lower dissipation. 
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Appendix A: Solution to Eq. (|55p 

In this Appendix we give some details on how to solve 
to Eq. ([55)) . We begin with the equation for Gf 

-ziJ+Gf-KzGf7^_ + i-G«Gf-^GfG«= ^^^^ 
-ievA^(G«f, -f,G^). 

Expressing H± in Eq. (jAip in terms of G^ by means of 
Eq. P^ and making use of the retarded component of 
the orthogonality condition ([SS]) 



(A2) 



Gl Gf + Gf G'^ = 



we bring Eq. (JA1[) to the form 

-zevA^(G«f,-f,G^O 
with the obvious notation [of. Eq. (|20l) ] 



(A3) 



(A4) 



Multiplying both sides of Eq. (jA3P by G;^ and making 
use of the normalization condition (G^)^ = 1, we arrive 
at 



Gf = 



evA^ 



^^+e+^/r 



(f, -G^f.G^). 



(A5) 



Consider now the equation for Gf 

1 



^iJ+Gf + iGf i7_ + ^(G^Gf - Gf G:^) (A6) 



2t 



1 



hB.r<K\ 



= -^evA^(G^f, - f,G^) - — (G^^Gf - GfG 

and the corresponding orthogonality condition 

G^ Gf + G'f^Gt + Gf G^ + Gf G:^ = 0. (A7) 

It is convenient to separate Gf into a "regular" part, 
similar in form to Eq. psp . and an "anomalous" one as 
follows: 

G^^Gl + Gl, G[^GfF_-F+G^, (A8) 

where F± = F(e±aj/2). A substitution of this definition 
into the orthogonality condition (|A7|) yields the much 
simpler equation 



a /AA 



Gf GJ + G? G^ = 



(A9) 



Substituting Eq. (|X8)) into Eq. ((M|) and using Eq. (PTj) 
we find the equation for Gj 



The solution to this equation can be found by following 
the procedure similar to that used above to solve for Gf . 
First, we express H+ and H- in terms of G^ and G^, 
respectively: 



-^^^ + ^)G^Gt-(-^(- + ^]GlG^ 



G'f {F+ f,-f,F^)- (F+ f, - f, i^_) G^ 



(All) 



-lev A,, 



Then, using the orthogonality condition Eq. (|A9[) . we 
eliminate G^ in favor of G^: 



{~ii+ - i£.^ + l/r)Gf G1 = -levA^ (A12) 

X Gf (/'+ f, - f, F_) - (/'+ f, - f, F_) G^T 



The solution of this equation is obtained by multiply- 
ing both sides by Gf and employing the normalization 
condition (Gf )^ = 1: 



Gl 



evA„ 



^f + e + z/r 



(A13) 



F^ f, -f,F_- Gf iF+ T, - f, F_)Gi 



Equations (|X5|) and (|A13j) are valid for Ej < A, 1/r 
and arbitrary relation between A and 1/t. In the dirty 
limit At ^ 1, Eqs. (|A5[) and (|A13[) can be respectively 
approximated as 

Gf = -ievA^T(f,-Gff,G^), (A14) 

G1 = -ievA^T (A15) 

X [f+t, -f,i^_-Gf (F+f,-f,^_)G^ . 



iH+Gl + iG1H^ + —{Gf Gl - Gl G^l) 

2t 



— ievA^ 



Gf (F+ f, - f, F_) - (i^+ f, - % F-)Gi 



(AlO) 



Substituting these equations into Eq. 
Eq. dSH). 
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